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1.0 introduction 


This document provides the level B/C mathematical specifications 
for the Area Targets and Space Volumes Processor (ATSVP) . Pur- 
suant to the requirements of reference 1, this processor is de- 
signed to compute the acquisition-of~signal (AOS) and loss-of- 
signal (LOS) times for the following: 

a. Area targets 

(1) Earth- referenced circles which are specified by a 
latitude, longitude, altitude, and radius. 

(2) Celestial circles which are specified by a right 
ascension, declination, and angular radius. 

(3) Earth- referenced polygons which are an arbitrary 
Earth- fixed figure having up to five sides with the 
"corner points" defined by latitude, longitude, and 
altitude. 

(4) Celestial polygons which are an arbitrary, inertially 
fixed figure having up to five sides with the corner 
points defined by right ascension and declination on 
the celestial sphere. 

b. Space volumes 

(1) Earth- referenced space volumes which are an arbitrary, 
Earth- fixed, five- sided polyhedron. These volumes are 
defined by a lower- limit polygon at an altitude, h^, 
and the projection of this polygon to an altitude, . 
The corner points of the polygon are defined by 
latitudes and longitudes and rotate with the Earth. . . 

(2) Celestial-fixed space volumes which are an arbitrary, 
inertially fixed, five-sided polyhedron. These volumes 
are defined by a lower- limit polygon at an altitude, h^, 
and the projection of this polygon to an altitude, h^. 
The corner points of the -polygon are defined by right 
ascension and declination on the celestial sphere. 
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The AOS and LOS times for these targets are defined (ref. 1) as 
follows : 

a. Ground circles and polygons 

AOS - the time corresponding to the first subsatellite point 
to lie just inside the area. 

LOS - the time corresponding to the last subsatellite point 
just prior to exiting the area. 

b. Celestial circles and polygons 

AOS - the time corresponding to the first zenith point to lie 
just inside the area. 

LOS - the time corresponding to the last zenith point just 
prior to exiting the area. 

c. Earth- referenced and celestial-fixed space volumes 

AOS - the time at which the spacecraft (S/C) is just entering 
the volume. 

LOS - the time just prior to the S/C exiting the volume. 

Six data tables will contain the information necessary to completely 
describe the area targets and space volumes. These tables (ref. 1) 
are as follows : 

a. Ground targets table containing 10 targets in 1 block of data. 

b. Celestial circles table containing 10 targets in 1 block of data. 

c. Ground polygons table containing 20 targets in 2 blocks of data. 

d. Celestial polygons table containing 10 targets in 1 block of data 

e. Earth- referenced space volumes table containing 10 targets in 
1 block of data. 

f. Celestial-fixed volumes table containing 10 targets in 1 block 
of data. 

Section 2 of this document describes the characteristics of the 
area targets and space volumes and provides the mathematical 
equations necessary to determine whether the S/C lies within 
the area target or space volume. These equations provide a 
detailed model of the target geometry and will be used during 
the precise numerical search. 
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Section 3 discusses a semianalytical technique for predicting 
the AOS and LOS time periods. This technique is designed to 
bound the actual visibility period using a simplified target 
geometry model and unperturbed orbital motion. Its principal 
purpose is to reduce the burden on the precise numerical search 
by eliminating regions of the S/C orbit where AOS and LOS times 
are physically impossible. Section 4 provides a functional over- 
view of the ATSVP. This section outlines the overall process 
required to determine precise AOS and LOS times. 

Section 5 presents the detailed logic flow for the ATSVP. This 
section integrates the functional overview presented in section 
4 with the equations and approach presented in sections 2 and 3 
and the appendixes. Appendix A discusses the procedure for sub- 
dividing complex concave polygons into two or more simpler convex 
segments. The purpose of this subdivision process is to permit 
the equations in section 2 to be used on a segment-by- segment 
basis to test for containment. Appendix B .provides a solution 
to the conic intersection equations used in section 3 for 
belestial-f ixed targets. 
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2.0 AREA TARGETS AND SPACE VOLUMES CHARACTERISTICS AND CONTAINMENT 
CRITERIA 

The following subsections discuss the characteristics of each of 
the area targets and space volumes presented in section 1. The 
mathematical equations necessary to determine whether the S/C lies 
within the area target or space volume are also developed and dis- 
cussed. 

Two reference coordinate systems will be used. The inertial 
Aries mean-of-1950 (M50) coordinate system (fig. 2-1) will be the 
reference system when dealing with area targets and space volumes 
which remain inertially fixed. The rotating geocentric coordinate 
system (fig. 2-2) will be used when dealing with area targets and 
space volumes which rotate with the Earth. 

Each of the following six subsections is further subdivided into 
three topics: 

a. Procedure - a brief description of the steps to be performed. 

b. Equations - a statement of the input parameter requirements 
and development of the mathematical equations. 

c. Assumptions and limitations - a description of any simplifying 
assumptions and/or mathematical restrictions. 

For convenience, section 2.7 summarizes the equations for all of 
the area targets and space volumes. 
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"NAME : 

.ORIGIN: 

ORIENTATION: 


CHARACTERISTICS: 


Aries mean-of-1950, Cartesian, coordinate .system. 

The center of the Earth. 

The epoch is the beginning of Besselian year 1950 or Julian 
ephemeris date 2433282.423357. 

The ■ x M _y M plane is .the mean Earth's equator of epoch. 

The X.. axis is directed towards the mean vernal equinox of epoch 
M 

The Zj^ axis is directed along the Earth's mean rotational axis 
of epoch and is positive north. 

The Y m axis completes a right-handed system. 

Inertial, right-handed, Cartesian system. 


Figure 2-1.- Aries mean-of-1950 coordinate system. 
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NAME: 

ORIGIN: 

ORIENTATION: 


CHARACTERISTICS: 


Geocentric coordinate system 
Center of the Earth 

X„ - Y_ plane is the Earth’s true-of-date equator 
G G 

X„ passes through the Greenwich meridian 

is along the Earth' s rotational axis 

Y_ completes the right-handed system 
G 

Rotating, right-handed. Earth-fixed 


Figure 2-2.- Rotating geocentric coordinate system. 



2.1 EARTH- REFERENCED CIRCLES 


Earth- referenced circles are defined to be circular ground target 
areas whose centers are defined by geodetic latitude, longitude, 
and altitude (fig. 2—3) . The S/C lies within this ground target 
area if its subsatellite point lies within the perimeter of the 
circular area. 

2.1.1 Procedure 

The following procedure will be used to determine whether the S/'C 
lies within the ground target area: 

a. The geodetic coordinates of the ground target area will be 
transformed to the geocentric system, 
b- The S/C position vector will be transformed from the M50 
system to the geocentric system, 

c. A test will be performed to determine whether the S/C subsat- 
ellite point lies within the perimeter of the ground target 
area. 

2.1.2 Equations 

The following parameters are required: 

4> and X - geodetic latitude and longitude, respectively, of the 

center of the Earth-referenced circle (fig. 2-4) 

h - altitude of the Earth- ref erenced circle, measured with 

respect to the Fischer ellipsoid of 1960 (fig. 2-4) 

r - radius of the Earth-referenced circle (fig. 2-3)' 
c 

it - S/C position vector in the M50 system (fig. 2-3) 

t - time corresponding to R gc 
rnr*T 

[RNPj^g -Rotation, nutation, and precession (RNP) matrix which 

is used to transform vectors from the M50 system to the 
true-of- epoch inertial (TEI) coordinate system 
t - epoch time corresponding to the RNP matrix 

R m “ mean equatorial radius for the Fischer ellipsoid of 1960 
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NAME : Geodetic coordinate system. 

ORIGIN: This system consists of a set of parameters rather than 

a coordinate system; therefore, no origin is specified. 

ORIENTATION: This system of parameters is based on an ellipsoidal model 

of the Earth (e.g. , the Fischer ellipse of 1960). For 
any point of interest we define a line, known as the geodetic 
local, vertical, which is perpendicular to the ellipsoid 
and which contains the point of interest. 

h, geodetic altitude, is the distance from the point of 
interest to the reference ellipsoid, measured along the 
geodetic local vertical, and is positive for points out- 
side the ellipsoid. 

X is the longitude measured in the plane of the Earth's 
true equator from the prime (Greenwich) meridian to the 
local meridian, measured positive eastward. 

4> is the geodetic latitude, measured in the plane of the 
local meridian from the Earth's true equator to the geodetic 
local vertical, measured positive north from the equator. 

NOTE: A detailed explanation of declination, geodetic 

latitude, and geocentric latitude is provided 
on figure 2-4 (b). 

CHARACTERISTICS: Rotating polar coordinate parameters. Only position vectors 

are expressed in this coordinate system. Velocity vectors 
should be expressed in the Aries mean-of-1950, or the Aries 
true-of-date, polar for inertial or quasi-inertial repre- 
sentations, respectively. The Fischer ellipsoid model 
should be used with this system. 

(a) Basic definitions. 

Figure 2-4.- Geodetic coordinate system. 
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NAME: Geodetic coordinate system of point P. 

DEFINITIONS: h is the altitude of point P measured perpendicular 

from the surface of the referenced ellipsoid. 

$ is the geodetic latitude of point P. 

4> c is the geocentric latitude of point P. 

6 is the angle between radius yector and equatorial plane 
(declination) . 

X is the longitude of point P. Angle (+ east) between 
plane of the figure and the plane formed by the Greenwich 
meridian. 


(b) Detailed explanation! 


Figure 2-4.- Concluded. 
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F - flattening coefficient for the Fischer ellipsoid of 1960 

w - Earth rotation rate 
e 


The first step is to transform the ground target area from the geo- 
detic system to the geocentric system. The vector from the 
center of the Earth to the center of the ground target area can 
be expressed in the rotating geocentric system by 


where 



1 

(h + a p ) cos X cos <{> 

' (h + a®) sin X cos $ 

[h + (1 - E) 2 a p ], sin <j> 


em 

Vcos 2 <f> +~ (1 - E) 2 sin 2 ^ / 


( 2 - 1 ) 


( 2 - 2 ) 


The second step is to transform the S/C position vector from "the 
M50 system to the- geocentric system. This is accomplished by 

a. Transforming the vector from the M50 system to the TEI system 
using the RNF matrix. 

b. Transforming the resultant vector from the TEI system to the 
geocentric system. 


The S/C position vector in the TEI system is given by 


t^ 1 = [RNP]™^ S (2-3) 

sc *- -*M50 sc 

The S/C position vector in the geocentric system (fig. 2-5) is 
given by 
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Figure 2-5.- Relationship between TEI and- 
rotating geocentric systems,. 
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sin AA 
cos AA 
0 


0 


(2-4) 


where 


R G 

sc 


cos AA 
-sin AA 
0 


0 

1 


G 

£TEI 

X\ 

SC 

TEI 


AA = to (t 
e 


V 


(2-5) 


The final step is to determine whether the S/C subsatellite point 

lies within the ground target area (fig. 2-3) . The angle from 
^■G 

Cg to the perimeter of the circular ground target area, y^, 
is given (in degrees) by 



180 

7T 


( 2 - 6 ) 


The angle from 



i to the S/C , 


Y s , | is given by 


Y S 


cos 


-1 


S G . $ G 

sc B 


sc 


5 G 


(2-7) 


The S/C lies within the ground target area if 


^S 



( 2 - 8 ) 


2.1.3 Assumptions and Limitations 

The following assumptions are implicit in the equations presented 
in section 2.1.2: 

a. The S/C geocentric subsatellite point is used to compute entry 
into the ground target area. 

b. The effects of polar nutation and precession from time t to 

0 

time t can be neglected. 
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is a 


c. The radius of the Earth- referenced circle, r , 

c' 

segment of arc. 
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2.2 CELESTIAL CIRCLES 


Celestial circles are defined to be circular areas on the celestial 
sphere (fig. 2-6). The center of this area target is defined by 
right ascension and declination. The criterion for a S/C to lie 
within this area is for the S/C zenith point to lie within the 
perimeter of the celestial circle. 

2.2.1 Procedure 

The following procedure will be used to determine whether the S/C 
lies within the celestial circle: 

a. The unit vector along the centerline of the celestial circle 
will be computed in the M50 system. 

b. The dot product between this vector and the S/C position 
vector will be formed to determine whether the S/C lies 
within the celestial circle. 

2.2.2 Equations 

The following parameters are required: 

a A and - the right ascension and declination, respectively, of 
the centerline of the celestial circle expressed in 
the M50 system (fig. 2-6) . 

- the celestial circle angular radius (fig. 2-6) . 

- S/C position vector in the M50 system. 

The unit vector from the center of the Earth to the center of the 

/\ 

celestial circle, C g , is given by 


1 

f COS 

a A 

cos 

S A 

A. 1 





CB = | 

sin 

a A 

cos 

S A 



sin 

5 A 
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Centerline of 
celestial circle 



2-13 






The angle between this vector and the s/C, y , is 

fa 


Y S 


cos 


-1 



(2-10) 


The S/C lies within the celestial circle if 


Ye < Yl 


( 2 - 11 ) 


2.2.3 Assumptions and Limitations 
None. 
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2.3 EARTH- REFERENCED POLYGONS 

The Earth-referenced polygon is defined to be an arbitrary planar 
figure having up to five sides (fig. 2-7) . This figure is fixed 
with respect to the rotating Earth. The corner points (i.e., 
vertices) of this polygon are defined by geodetic latitude, longi- 
tude, and altitude. The basic criterion for penetration is to 
ensure that the S/C subsatellite point lies within the perimeter 
of the ground target area. 


2.3.1 Procedure 

The basic procedure for this ground target area is similar to the 

procedure presented in section 2.1.1. It consists of 

a. Transforming the geodetic coordinates of each polygon vertex 
to the geocentric system. 

b. Transforming the S/C position vector from the M50 system to 
the geocentric system. 

c. Testing to determine whether the S/C is interior to all 
planes defined by the sides of the polygon. 


2.3.2 Equations 

The following parameters are required: 


n 

d> and X . 

T i l 


h i 


R 


sc 

t 


[RNP] 


TEI 

M50 


em 


R 

€ 

F 


number of sides 

geodetic latitude and longitude, respectively, of 
each vertex (fig. 2-7) 

geodetic altitude of each vertex (fig. 2-7) 

S/C position vector in the M50 system 

» ■ — V" 

time corresponding to R 

s c 

RNP matrix to transform from the M50 system to the 
TEI system 

epoch time corresponding to the RNP matrix 
mean equatorial radius for the Fischer ellipsoid of 1960 
flattening coefficient for the Fischer ellipsoid of 1960 
Earth rotation rate 
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Greenwich 

meridian 


Figure 2-7.- Earth-referenced polygon. 
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The first step is to transform the parameters defining each 
vertex of the polygon from the geodetic system to the rotating 
geocentric system. These vectors are given by 

(im + a^) cos X ^ cos 
<h.j_ + a F ) sin cos 
[h ± + (1 - F) 2 aj sin 

where 

is defined by equation 2-2 with replacing <j>. 

The S/C position vector in the geocentric system, , is then 

sc 

obtained by using equations 2-3 through 2-5. 

The next step in the procedure is to determine whether the S/C 
lies interior to all planes defined by the sides of the polygon. 
Figures 2-8 and 2-9 illustrate the geometry. All vectors in these 
figures are with respect to the geocentric system. The centroid 
of the ground target area, C G , (fig. 2-8) is defined as 





(2-13) 


The unit normal vectors to each side of the polygon (fig. 2-9) are 
given by 



.n-1 (2-14a) 



(2-14b) 
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Figure 2 8.- Centroid of Earth-referenced polygons 
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the Earth 


Figure 2-9.- Containment test for Earth-referenced polygons. 
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(2-15) 


v Q 

The perpendicular distance from C B to each polygon side is 


a i - - 2 b) 


i • ill 


The S/C must lie interior to the x plane if 


"i • ?B i *i 


(2-16) 


where 


-*-G -±G _ ^G 


p„ = R 

sc 


B 


(2-17) 


If equation 2-16 is satisfied for all sides, then the S/C sub- 
satellite point lies within the perimeter of the polygon. 

2.3.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 
equations presented in section 2.3.2: 

a. The S/C geocentric subsatellite point is used to determine 
entry into the ground target area. 

b. The effects. of polar nutation and precession , from time t 

© 

to t can be neglected. 

c. The vertices of the polygon are specified in a counterclockwise 
order as viewed from the top. 

d. The polygon is convex; i.e., the interior angles between the 
sides defining the vertices are less than 180 degrees. 1 

e. The angular separation between two consecutive vertices is 
sufficient to permit a nonzero cross product (eq. 2-14) . 

f. The S/C and the ground target area lie in the same hemisphere. 
The procedure discussed in section 3 will assure this 
condition. 

1 Concave polygons can be accommodated by subdividing them into 
two or more convex polygons. Appendix A discusses this procedure. 
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2.4 CELESTIAL POLYGONS 


Celestial polygons are defined to be arbitrary figures having up 
to five sides with the corner points (i.e., vertices) defined by 
right ascension and declination on the celestial sphere. The 
criterion for - penetration into this area is to ensure that the 
S/C zenith point lies within the confines of the polygon. Figure 
2-10 illustrates the celestial polygon. It is noted that this 
polygon also represents a variable area polyhedron which remains 
inertially fixed. 


2.4.1 Procedure 

The basic procedure for this area target is similar to the proce- 
dure presented in section 2.3.1, i.e.,: 

a. The unit vector along the "centroid" of the polyhedron will 
be computed in the M50 system. 

b. Tests will be made to determine whether the S/C is interior 
to all planes defined by the sides of the polyhedron. 

2.4.2. Equations 

The following parameters are required: 
n - number of sides 

a i and - right ascension and declination, respectively, of each 
vertex in the M50 system (i = l,2,3,...n) 

3^ sc - S/C position vector in the M50 system 


The first step is to determine the unit vector along the centroid 
of the polyhedron. The unit vectors from the center of the Earth 
to each vertex in the M50 system are given by 


cos a. cos 6. 

i l 

R-. = { sin a . cos 6 . 
i I x x 

sin 6 . 
x 


> i = 1,2,3, ...n (2-18) 
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The unit vector along the centroid of the polyhedron is 


C 



where 



(2-19) 


( 2 - 20 ) 


The final step is to determine whether the S/C lies interior to. all 
planes defined by the sides of the polyhedron. Figure 2-11 illu- 

A 

strates the geometry. The slant range from C to the S/C is 

P B = B sc - £ (2-2D 

Equations 2-14 through 2-16 (with R. , C^, and p 0 replaced by 

„ 1 13 13 

R^, C, and p , respectively) are then used to determine whether 
the S/C is interior to all polyhedron planes. 


2.4.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 

equations presented in section 2.4.2: 

a. The vertices of the polygon are specified in a counterclockwise 
order as viewed from the top. 

b. The polygon is convex; i.e., the interior angles between the 

. . . 1 
sides defining the vertices are less than 180 degrees. 

c. The angular separation between two consecutive vertices is 
sufficient to permit a nonzero cross product (eq. 2-14). 

d. The S/C and the area target lie within the same hemisphere. 

The procedure outlined in section 3 will assure this condition. 


Concave polygons can be accommodated by subdividing them 
into two or more convex polygons as discussed in appendix A. 
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Figure 2-11.- Containment test for celestial polygons. 


2.5 EARTH- REFERENCED SPACE VOLUMES 

The Earth-referenced space volumes are arbitrary five-sided poly- 
hedrons which are defined by a lower-limit polygon at an altitude, 
h 1 , and the projection of this polygon to an altitude, The 

vertices of this polygon are defined by geodetic latitude, and 
longitude and rotate with the Earth. Figure 2-12 illustrates 
this type of space volume. As shown, the planar cross sectional 
area of this polyhedron remains constant with respect to altitude. 

2.5.1 Procedure 

The following procedure will be used to determine whether the S/C 
lies within the space volume: 

a. The geodetic parameters defining each vertex of the lower 
boundary will be transformed to geocentric position vectors. 

b. The centroid vector to the lower boundary will be computed 
in the geocentric system. 

c. The S/C position vector will be transformed from the M50 
system to the geocentric system. 

d. Tests will be performed to ensure that the S/C lies above the 
lower boundary and below the upper boundary. If either of 
these tests fails, then no further computations are re- 
quired. 

e. Assuming the previous step is passed, tests will be performed 
to determine whether the S/C is interior to all planes 
defined by the sides of the polyhedron. 

2.5.2 Equations 

The following parameters are required: 

<j) and X i - geodetic latitude and longitude, respectively, of 
each vertex of the lower boundary (i = 1,2 ,3, 4, 5) 
h^ - geodetic altitude of the lower boundary 

h 2 - geodetic altitude of the upper boundary 
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Space volume centerline 



Figure 2-12.- Constant area polyhedron. 


2-2 6 




R 


sc 


[RNP] 


TEI 

M50 


R 


em 


.e 


S/C position vector in the M50 system 
time corresponding to R 

s c 

RNP matrix to transform from the M50 to TEI system 
epoch time corresponding to the RNP matrix 
mean equatorial radius for the Fischer ellipsoid 
of 1960 

flattening coefficient for the Fischer ellipsoid 
of 1960 

Earth rotation rate 


The first step is to transform the parameters defining each 
vertex of the lower boundary to geocentric position vectors. 
Equation 2-12 (with h^ = 0) provides the necessary transformation. 


The centroid of these vectors is given by 



*G _ i=l 
u n 


( 2 - 22 ) 


The geocentric vector to the centroid of the lower boundary 
(fig. 2-13) is given by 



C G + 



(2-23) 


Next, the S/C position vector in the geocentric system, , is 

computed using equations 2-3 through 2-5. 

The fourth step in the procedure is to ensure that the S/C lies 
between the upper and lower boundaries. Figure 2-13 illustrates 
the geometry. All vectors in this figure are with respect to 
the geocentric system. The slant range from C g to the S/C, p B , 
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Centerline 



Figure 2-13.- Boundary tests for constant area polyhedrons. 
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is given by 





(2-24) 


The S/C lies above the lower boundary and below the upper boundary 
if 


h 2 ” h l - P B 


®b G 

®*bI 


> 0 


(2-25) 


Assuming that equation 2-25 is satisfied, the final step is to 
determine whether the S/C lies interior to the side planes of 
the polyhedron. 'Figure 2-14 illustrates the geometry. The unit 

y g 

normal vectors from C ' to each side of the polygon are given by 

B 


n g 


2 b x 




I«B X 




*?+x)l 


i = 1,2, 3, 4 (2-26a) 


Nr- = 


->G / ->*G 

X R^ 
B \ 5 




,-K3 /+ G 

|Cg X K 


{n - n) 


(2-26b) 


Equations 2-15 and 2-16 are then used to determine whether the 
S/C is contained within the polyhedron. 


2.5.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 
equations presented in section 2.5.2: 

a. The altitude of the lower boundary, h , is measured with 
respect to the centroid of the vertices which lie on the 
Fischer ellipsoid. 
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b. The effects of polar nutation and precession from time t 
to t can be neglected. 

c. The vertices are specified in a counterclockwise order as 
viewed from the top. 

d. The planar area of the polyhedron is convex; i.e. , the inter- 
ior angles between the sides defining the vertices are less 
than 180- degrees. ^ 

e. The local horizon for determining whether the S/C is between 
the upper and lower boundaries is perpendicular to the centroid 
vector. 

f. The distance between two consecutive vertices is sufficient 
to permit a nonzero cross product (eq. 2-26) . 


■'"Concave polyhedrons can be 
into two or more convex portions 


accommodated by subdividing them 
as discussed in appendix A. 
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2.6 CELESTIAL-FIXED SPACE VOLUMES 

The celestial-fixed space volumes are arbitrary five-sided poly- 
hedrons -.which are defined by a lower- limit polygon, at an altitude,, 
h^, and a projection of this polygon to an altitude, h 2 -' The 
vertices of the polygon are defined by, right ascension and de- 
clination and are inertially fixed. Figure 2-12 -can also be- 
used to illustrate this type of space volume., As : mentioned! 
previously, the planar cross sectional area of this- polyhedron- 
remains! constant with respect to altitude 

2-.6..1 Procedure 

The basic procedure for this space volume is very similar to the- 
procedure presented in section 2.5.1. It consists of: 

a. Computing the centroid vector to the lower' boundary in the 
M50 system. 

b. Testing to ensure that the S/C lies above, the. lower boundary 
and below the- upper boundary.. If either of these tests fails, 
then no 'further computations are. required. 

c. Assuming the previous step is. passed',, further tests; will be 
performed to determine whether, -the S/C is interior to all 
planes defined by the sides of the polyhedron. 

2.6.2 Equations. 

The following parameters are required:- 

a. and 6^ - right ascension and declination,, respectively, of 

each vertex of the lower boundary in. the M50 system- 
(1 = 1,2, 3,4,5) 

h - altitude of the lower boundary 

h 2 - altitude of the upper boundary 

R - mean equatorial radius 

em 
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The unit vector to the centroid of the polyhedron, C, is com- 
puted via equations 2-18 through 2-20. The vector to the. centroid 
of the lower boundary is given by 

S B " ( R ern + h l) 2 < 2 ' 27) 

Equations 2-24 and 2-25 (with variables R and replaced by 

+ + sc B 

R sc and Cg, respectively) are then used to ensure that the S/C 
lies between the upper and lower boundaries. 

Assuming equation 2-25 is satisfied, equation 2-26 (with 
and R? replaced by C_, and R. , respectively) is used to 
define the unit normal vectors to each side of the polyhedron. 
Finally, equations 2-15 and 2-16 are used to determine whether 
the S/C lies within the space volume. 

2.6.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 
equations presented in section 2.6.2: 

a. The altitudes of the lower and upper boundaries are measured 
with respect to the mean equatorial radius. 

b. The vertices are specified in a counterclockwise order as 
viewed from the top. 

c. The planar area of the polyhedron is convex; i.e., the in- 
terior angles between the sides defining the vertices are 

n 

less than 180 degrees. 

d. The distance between two consecutive vertices is sufficient 
to permit a nonzero cross product (eq. 2-26) . 


Concave polyhedrons can be accommodated by subdividing 
them into two or more convex portions as discussed in appendix A. 
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2.7 SUMMARY OF EQUATIONS 

This section summarizes all of the equations presented in the 
previous sections for the various area targets and space volumes. 
The order of presentation and equation numbers correspond to the 
computation sequence discussed in the text. 


2.7.1 Earth- Reference d Circles 

(h + a cos A cos <f> 
^h + a^ sin A cos <p , 
[h + (1 - F) 2 a F ]sin <f> 

where 



(■ 2 - 1 ) 


a F 


V 


R 

em 

2 ; 2 2 ' 
cos <p + (1 - F) sin 


( 2 - 2 ) 


^TEI _ [gjjp-j 
SC L J 


TEI ± 
M50 sc 


(2-3) 


R 

sc 


coe AA 
-sin AA 
0 


sin AA 
cos AA 
0 


0 

0 

1 


G 


£TEI 

K. 

SC 


JTEI 


(2-4) 


where 



(2-5) 
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(2-6) 


Y 


A 



180 

7T 


Y S = 





2.7.2 Celestial Circles 


(2-7) 


( 2 - 8 ) 


/\ 



{ cos a cos 
sin a cos 



Y q = cos 


- 1 - 


sc 


"B 


i -*■ 

; R 


sc 


Y < Y, 
! S - t A 


2.7.3 Earth-Referenced Polygons 


(2-9) 


( 2 - 10 ) 


( 2 - 11 ) 



( ( h i + a F ) 

cos X . 

l 

COS 

*i) 



R? = 

l 

( h i + a F ) 

sin X i 

COS 

+i 

i = 1,2, 3, . . .n 

(2-12) 


([ h i. + (1 - 

5 

NJ 

P 

1 1 

sin 

+i) 
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where 


R 


a F 


era 


^cos 2 cp i + (1 - F) 2 sin 2 ^ 


R^“ =■ [RNPlr-Kn R 
sc J M50 sc 


->G 

R 

sc 


cos AX sin AX 0 
-sin AX cos AX 0 
0 0 1 


£TEI 

sc 

TEI 


where 


AX - <o a (t - t e ) 




t 

i=l 


n 


1 


l*? + i x *2' 


= 1,2,3,., .n-1 


N G = 
n 


S? X R G 


n 


i ->G -*G i 

|R£ x R~| 


d. 

l 


= »?•• [n - %) i - 1^-3,. 


. .n 


Ct G -*g „ , 
N. • < d. 

l K B — l 


i = 1,2,3, n 


( 2 - 2 ) 

(2-3) 

(2-4) 

(2-5) 

(2-13) 

(2-14a) 

(2-14b) 

(2-15) 

(2-16) 
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where 


n' G = r g - r G 
P B R sc C B 


(2-17) 


2.7.4 Celestial Polygons 


R i “ 


/ cos a . cos £ . 


sin a. cos 6. 


sin 6 . 

i 


i = 1,2,3,. ..n (2-18) 


where 


C = 


I C| 


(2-19) 


C = 


n 

_ A 

E R i 

i=l 


n 


(2-s20) 


= R - C 
B sc 


( 2 - 21 ) 


A R. X R. 

N. = 1+x — i = 1,2,3, ...n-1 (2-14a) 

1 R. .. X R. 

x+1 i 


N = 


R. X R 
1 n 


11 |E l XE nl 


(2-14b) 


di = N ± • (L - cj i = 1,2, 3, ...n (2-15) 


N 1 ' p B i d l 


i = 1,2,3, ...n (2-16) 
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2.7.5 Earth-Referenced Space Volumes 


where 




a F 

cos A^ 

cos 



S F 

sin A . 
r 

cos 

♦i 

► i = 1,2,3,4,5(2-12) 

(1 ' 

" F) 2 a F 

sin 

^ii 

| 


R 


a = 
F 


em 


cos 2 *, + (1 - F) 2 sin^tj) . ^ 
1 1 


( 2 - 2 ) 


E *i 

S G = ±=l 


(2-22) 


G B G n l i+G, 


(2-23) 



+TEI 

sc 

= 


-V 

R 

sc 


(2-3) 


cos 

AA 

sin AA 

0" 

G 


r g = 
sc 

-sin 

AA 

cos AA 

0 

r tei 

sc 

(2-4) 

0 


0 

1 . 

TEI 



where 


AA = o) 


(* - 0 


(2-5) 


o G = R G - c G 
P B SC C B 


(2-24) 
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(2-25) 


h 2 - h x > 


S B 

0 

I c bI 


S B X 

l J B X 


- S i + l) 


S9, 1 
1 + 1 


i = 1, 2,3,4 (2-26a) 


_ 2 B X (*5- S ?) 


2 b x (*< 


*S)I 


(2-26b) 


. -G /*G £ g \ 

1 . = N . • R . - 

l l \ x B/ 


i = 1,2, 3, 4, 5 (2-15) 


fi ? • ^B i d i 


i =•' 1,2, 3, 4, 5 (2-16) 


2.7.6 Celestial-Fixed Space Volumes 


cos a . cos 6 . 

l l 

A 

R. = sin a . cos 6 . 

i li 


sin 6 


i = 1,2, 3-4, 5 (2-18) 


(2-19) 


5 A 

E R i 

i=l 

5 


(2-20) 


C = (r + h ) C 
B \ em 11 


(2-27) 
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(2-24) 


0 (2-25) 

■ i = 1,2, 3, 4 (2-26a) 

(2-26b) 

i = 1,2, 3, 4, 5 (2-15) 
i = 1,2, 3, 4, 5 (2-16) 



3 . 0 SEMIANALYT ICAL ALGORITHM TO PREDICT AOS AND LOS TIMES 


This section presents a semianalytic technique for predicting AOS 
and LOS tiroes. This technique is designed to bound the actual 
visibility periods and uses a simplified geometry model in con- 
junction with unperturbed orbital motion. Its principal purpose 
is to reduce the burden on the precise numerical search by 
eliminating regions of the S/C orbit where AOS and LOS times are 
physically impossible. 

The basic procedure consists of the following: 

a. Determining the parameters of an outer boundary cone which 
circumscribes the entire area target or space volume. 

b. Determining whether the orientation of the S/C orbit plane 
will permit any intersections of the outer boundary cone 
with the S/C orbit plane. 

c. Assuming intersections are possible for the following: 

(1) Celestial- fixed targets 

(a) Compute the conic intersection points of the outer 
boundary cone with the S/C orbit plane. 

(b) Compute predictions of AOS and LOS times 

Ct___ and t rne , , respectively) for each 

revolution based upon these intersection points. 

{2) Earth- fixed targets 

(a) Compute the time of closest approach (TCA) of the 
S/C position vector with the centerline of the 
outer boundary cone and the corresponding closest 
approach angle (y^) . These computations are 
performed on a revolution-by- revolution basis. 

(b) If is less than or equal to the half-cone 

angle of the outer boundary cone, then compute 
predictions of AOS and LOS times based upon the 
times when the S/C lies on the perimeter of the 
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outer boundary cone. These times are denoted by 
t AOS and t LOS * 

jS. j\ 


After the AOS and LOS time predictions are determined, the precise 
numerical search (using the detailed area target and space 
volume geometry ■ model discussed in sec. 2) can then be performed 
over the reduced time regions from t 


AOS 


to 


K 


'LOS 


K 


The following parameters are -required to predict the AOS and LOS times 


t start 

^nd 

At 



[R3SJp] TE ^ 

lj M50 


t 

e 


- start time for the search 

- end time for the search 

- time increment for the search. For the purposes of 
this algorithm. At will be used as a tolerance for 
converging upon the AOS and LOS times for Earth- 
'referenced targets. 

- S/C position vector in the M50 system at t . . . 


- S/C velocity vector in the M50 system at t start 

- RNP matrix to transform from the M50 to the TEI 
system. 

- epoch time corresponding to the RNP matrix. 


The following subsections provide the necessary equations. 

Section 3.1 discusses the computation of the outer boundary cone 
parameters. Section 3.2 discusses the technique to determine 
whether any AOS or LOS times are possible. Section 3.3 presents 
the equations to predict AOS and LOS times. Section 3.4 discusses 
the assumptions and limitations implicit in these equations. 

Two reference coordinate systems are used in these computations. 
The M50 coordinate system is used when dealing with celestial- 
fixed targets and the TEI system is used for Earth- referenced 
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targets. The TEI system was selected for Earth-referenced targets 
in order to model the effects of polar nutation" and precession 
on the Earth's spin axis. 


For celestial-fixed targets, the M50 Keplerian elements 

i f fi, w and M ) at the start time will be computed from 
° 

and V 

sc o 



e,, 



For Earth-referenced- targets, the TEI Keplerian elements at the 

-+TET 

start time will be computed from R .and . Where 

Sc o 


+TEI 

sc 

o 


s !! 1 = [*»:£? * 


sc 


M50 sc 


(3-1) 


- CrnpD^ 


sc 


J .M50 sc 


(3-2) 
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3.1 COMPUTING OUTER BOUNDARY CONE PARAMETERS 

The outer boundary cone is defined by 

a - right ascension of the centerline 
o 

<b - latitude of the centerline 
T o 

- half-cone angle 

For Earth-referenced targets these parameters will be defined in 
the TEI system at t t"* For ce l es ' t: ^ a l -: ^ ; *- xe< ^ targets, these 

parameters will be defined in the M50 system and remain invariant 
with respect to time. 

3.1.1 Earth- Referenced Circles 

The unit vector along the centerline of the cone is given in the 
TEI system by 

cos AA - sin AA 0 ’ TEI -+G 
C TEI = sin AA cos AA 0 -A_ (3-3) 

L o o iJ G I^B I 

is given by equation 2-1 

is given by equation 2-5 with t equal to t 

S"t3.3T*C 

Thus, the right ascension and latitude of the centerline is given 
in the TEI system (at the start of the search) by 

pTEI- 

zz 

C TEI 
x 

^The TEI system is selected for Earth-referenced targets in 
order to model the polar nutation and precession from M50 reference 
to the epoch time of the RNP matrix. 


(3-4) 
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where 


"TEI . / TEI TEI TEl\ 

C ' \ x ' y ' L z / 

The half-cone angle, y a * is given by equation 2-6. 

3.1.2 Celestial Circles 

The right ascension, latitude, and half-cone angle in the M50 
system are given via input 


“o = “a 

(3-6) 

*o = S A 

(3-7) 


3.1.3 Earth-Referenced Polygons . 

The unit vector along the centerline of the outer boundary cone is 
given in the TEI system by equation 3-3 where equation 2-13 is used 

-4-r* 

to compute C„. Equations 3-4 and 3-5 are then used to compute 

Jd 

.a and <f> 0 , respectively. The outer boundary cone is centered 
along the centroid vector and has an angular radius which contains 
the vertex furthest from the centroid vector (fig. 3-1). Thus, 
the half-cone angle is given by 

(3-8) 


Y a = max 


cos 


-1 


n • 

% 

|R?I 

iagi 


i = 1,2,3, . . .n 


max implies the maximum algebraic value 
is given by equation 2-12. 


where 



Cone 

circumscribing 

Earth-referenced 

polygon 


Centerline of 



Figure 3-1.- Circumscribing polygons. 
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3.1.4 Celestial Polygons 

The unit vector along the centerline of the celestial polygon is 

v A ipg J 

given by equation 2-19. - Equations 3-4 and 3-5 (with C replaced 
by C) are then used to compute a Q and <j> o ; respectively. The 
half-cone angle which contains the vertex furthest from the centroid 
is given by 



i = 1,2,3, 



• (3-9) 


where 

A 

R/ is given by equation 2-18. 


3.1.5 Earth-Referenced Space Volumes 

The unit vector along the centerline of the' outer boundary cone 
for Earth- referenced space volumes is given by equation 3-3 with 

Cg replaced by C (eq . 2-22). Equations 3-*4 and 3-5 


are then used to compute a o and <f> 0 ' respectively. Since the 

Earth-referenced space volume represents a constant area polyhedron 
(sec. 2.5), the outer boundary cone to be used for predicting AOS 
and LOS times circumscribes the space volume planar area at the 
S/C altitude. However, for noncircular orbits, the S/C altitude 
varies as a function of time and thus the size of the outer 
boundary cone would also vary as a function of time. The size of 
the cone is largest at perigee and smallest at apogee. Since the 
purpose of the outer boundary cone is to' bound the actual visibility 
region, the half-cone angle at the S/C perigee altitude will be 
computed and held constant with respect to time. The vector from 


the center of the Earth to each vertex of the space volume pro- 


jected to the S/C perigee altitude, R^ , (fig. 3-2) is given by 



i = 1,2,3, . . .n (3-10) 
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Figure 3-2.- Circumscribing space volumes. 
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where 


± Q 

R7 . is given by equation 2-12 

"^G 

C is given by equation 2-23 
£ 

R is the S/C perigee radius 

Xr 

- a(l - e) (3-11) 

Thus the half-cone angle which contains the vertex furthest from 
the centerline of the cone is given by 


Y , = max 





cos -1 

1 B 

i 1,2,3,. » . n 


l*?l 1^1 

. 


(3-12) 


3.1.6 Celestial-Fixed Space Volumes 

The unit vector along the centerline of the outer boundary cone 
for celestial- fixed space volumes is given by equation 2-19. 

Arrrn j a 

Equations 3-4 and 3-5 (with C replaced by C) are then used 
to compute a Q and <|> 0 , respectively. Since the celestial-fixed 

space volume also represents a constant area polyhedron (sec. 2.6), 
the outer boundary cone to be used for predicting AOS and LOS times 
circumscribes the space volume planar area at the S/C perigee 
altitude (sec. 3.1.5). The vector from the center of the Earth to 
each vertex of the space volume at the S/C perigee altitude, 

-y ' . . 

R^ , is given by 


-V 



R. + 
1 




(3-13) 
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where 


Thus 


A 

is given by equation 2-18 

2 is given by equation 2-27 
B 


r is the S/C perigee radius 

p 

(equation 3-11) 




max 


-1 


cos 



c 


i 1^2f3f«* *n 


(3-14) 
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3.2 DETERMINING WHETHER AOS AND LOS TIMES ARE POSSIBLE 


The following computations will determine whether the relative 
orientation of the S/C orbit plane and the outer boundary cone 
will permit any intersections between these two figures. If no 
intersections are possible, then no AOS or LOS times are possible. 
Since the Earth-referenced targets rotate with the Earth and the 
celestial-fixed targets remain inertially fixed, two' different 
types of tests are required. The celestial- fixed situation is 
the simplest and is discussed first. 

3.2.1 Celestial-Fixed Targets 

Figure 3-3 illustrates a celestial-fixed outer boundary cone and 
a S/C orbit plane.- The S/C orbit plane is perpendicular to the 
unit angular momentum vector which, is given by 


A 


h = 



(3-15) 


where 

R and ^ are the S/C position and velocity vectors, re- 
SC sc 

u spectively, in the M50 system at the beginning 

of the search. 


The centerline of the outer boundary cone, C, is given in the 
M50 system by 


\ 


cos 

a 

cos <j> 

o 

o 

sin 

a 

cos d> 

o 

r o 


sin J 


(3-16) 


3-11 




plane 


Figure 3-3.- Intersection of celestial-fixed 
outer boundary cone with S/C orbit plane. 
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where 


a Q and. <j> Q were computed using the procedure discussed in 
section 3.1. 

The intersection points of the outer boundary cone with the S/C 

A A 

orbit plane are denoted by 1^ and on figure 3-3. These 

vectors can be found by solving the following conic intersection 
equations'^ 

h • I — 0 (3-17a) 

c • I = cos y a {3-17b) 

1 . 1=1 (3-17c) 

A 

These equations provide 0/ 1/ or 2 solutions for I. If no 
solutions are produced, then the relative orientation of the 
outer boundary cone and the S/C orbit plane does not permit 
intersections, and no AGS or LOS times are possible. If only 
one solution is produced, the AOS and LOS times are the same. 

If two solutions are produced, the AOS and LOS times predictions 
are made using the method outlined in section 3.3, 


The method for solving these equations is presented in 
appendix B. 
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3.2.2 Earth- Referenced Targets 


Figure 3-4 illustrates an Earth-referenced outer boundary cone 
at a particular instant of time. This cone rotates about the. 
Earth's polar axis. Hence, over one sidereal day, the voltime of 
space swept out by the outer boundary cone describes a spherical 
sector whose upper. and lower limits are determined by y and 
<j> o (fig. 3-5) . Intersections of the S/C orbit plane with .this 
volume of space are possible only if the S/C orbit inclination' 
exceeds the lower limit of the spherical sector. 

For posigrade and polar orbits, the S/G orbit plane will intersect 
the spherical sector only if 



(3-18a) 


For retrograde orbits, intersections are- possible only if 


1 


180 - i-> |# 0 | - y K 


(3-18b) 


where 

i is the S/C orbit inclination in the TEI system 

If equation 3-18 is not satisfied, then no intersections' are 
possible 1 and thus no AOS or LOS times are possible. If equation 
3—18 is satisfied, then the approach presented i-n section 3.3 is 
used to predict AOS. and LOS times. 


The absolute . value sign is used' in these equations to ac- 
commodate targets in both the northern and southern hemispheres. 
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Figure 3-4 . - Earth- ref erenced outer boundary cone . 



Upper limit 



Figure 


3-5.- Spherical sector swept out by Earth 
referenced outer boundary cone. 
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3.3 PREDICTING AOS AND LOS TIMES 

The following subsections present the equations necessary to 
predict AOS and LOS times for celestial-fixed and Earth-referenced 
targets. These equations assume that the tests performed in 
section 3.2 have indicated that intersections of the outer boundary 
cone with the S/C orbit plane are possible. 

3.3.1 Celestial-Fixed Targets 

/V A 

The intersection vectors, 1^ and I (sec. 3.2.1), can be 
used directly to predict the AOS and LOS times. Figure 3-6 
illustrates the geometry. The unit vector in the direction of 

A 

the ascending node, ft, is given by 


ft 


cos ft 
sin ft. 


(3-19) 


where 

ft = right ascension of the ascending node in the M50 
system 

A 

The angles between ft and the two intersection, vectors are given 
by 


u 


1 


tan 


h • (ft X I ) 

A A 

ft * I 1 


U£ = tan 


-1 


[h • (ft X l 2 ) 


/V A 

ft * I. 


(3-20a) 


(3-20b) 
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z 


M 



Figure 3-6.- Definition of intersection points. 


3-18 


One of these angles corresponds to the AOS point and the other to 
the LOS point. The appropriate angle can be determined by noting 
that the visibility period must be less than or equal to one half 
of the orbit. 


Thus let 


K = tan 


-1 


h . (I 1 X l 2 ) 

A A 

I 1 * X 2 


(3-21) 


If K > 0 


If K < 0 


The true anomalies of 


U AOS 

= u. 

U LOS 

= u, 

U AOS 

= u, 

U LOS 

— u. 


the AOS and LOS points 


( 3- 22a) 
(3-23a) 

(3-22b) 

(3-23b) 

are given by 


f aos/los u aos/los ” 03 


(3-24) 


where 

(d is the argument of perigee of the S/C orbit. 
These can be converted to AOS and LOS times by 


'AOS/LOS 


_ m aos/los 


+ T 


(3-25) 


M 


where 


m aos/los e aos/los e sin e aos/los 


(3-26) 
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e aos/los 2 tan 


-1 


tV 


1 - e' 
1 + e 


tan f - AOS/MB 


( 3 - 27 ) 


M = 


360 


(in degrees per unit of time) 

T - 2 *Jv- 


t = t 


start 


M . 

■* ~ + (R - 1)- T 
M 


( 3 - 28 ) 

( 3 - 29 ) 

( 3 - 30 ) 


R is the revolution number measured with respect to 
^start 

a is the semimajor- axis of the S/C orbit 
e is the eccentricity of the S/C orbit 
M q is the S/C mean anomaly at ^gtart 
p is the gravitational constant 


Since both the orbit plane and the celestial-fixed target remain 
inertially fixed, u AQS and u LQS remain constant. Thus, the 
AOS and LOS times for future revolutions differ by only the 
orbital period. The AOS and LOS times for- the subsequent N 
revolutions are given by 

^"AOS/LOS K+1 = WlOS k + T K “ 1 ' 2 ' 3 "‘ N (3 - 31) 
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3.3.2 Earth-Referenced Targets 

Computation of AOS and LOS times for Earth-referenced targets is 
somewhat more complicated than for celestial-fixed targets since 
the orientation of the Earth-referenced target varies as a func- 
tion of time. The following procedure will be used to estimate 
the AOS and LOS times : 

a. The TCA of the center of the Earth-referenced outer boundary 
cone with the S/C position vector will be determined along 
with the corresponding closest approach angle, y CA . These 
computations will be performed on a revolution-by-revolution 
basis . 

b. The closest approach angle from the previous step will be N 
compared with y A * if Y CA > y a , then no AOS or LOS times 
are possible for that revolution, and the previous step will 
be repeated for the next revolution. If y^ A <_ y , then 
AOS and LOS times will be computed by solving for the times 
when the angle between the center of the outer boundary cone 
and the S/C position- vector is equal to y . 

3 . 3 . 2 . 1 Determining TCA 

The closest approach point occurs when the S/C position vector 
: and the center of the outer boundary cone lie in the same hemis- 
phere and are coplanar (fig. 3-7). Both the S/C and the center 
of the outer boundary cone are moving with time. The center of 
the cone is rotating about the Earth's Z axis at the Earth's 
rotation rate. The S/C is moving about the orbit angular momen- 
tum vector at a rate which is dependent upon both the orbit period 
and eccentricity. As a result, the TCA cannot be analytically 
determined for noncircular/nonequator ial orbits. The TCA can, 
however, be determined iteratively by the following: 
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Figure 3-7.- Closest approach point. 
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a. Estimating the closest approach point (for the first revolu- 
tion, this estimate would be based upon t . . ; for sub- 

s tdrt 

sequent revolutions, this estimate would be based upon the 
TCA of the previous revolution plus one orbital period) . 

b. Determining the time at which the S/C will reach this point. 

c. Updating the closest approach point based upon the time 
estimate from the previous step. 

d. Repeating the two previous steps until the difference in 
time between two successive iterations is less than or equal 
to the time step. At. 


The following discussion develops the equations necessary to 
compute the TCA. 

The unit vector from the center of the Earth to the center of 
the outer boundary cone is given by 


/ . \ 


- 1 

cos 

a 

cos 


C = | 

f sin 

a 

cos 

to 


sm 

<j> 

Y o 



(3-32) 


where 


a — a + 
o 


“e (t - t start ) 


(3-33) 


a Q is the initial right ascension of the centerline 
(sec. 3.1) 


o> e is the Earth rotation rate 

<j> Q is the latitude of the centerline (sec. 3.1) 
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The closest approach point, r^, occurs when 


A A 

: _ h x c 

‘ca ~ t r 77 

h x C 


x h 


(3-34) 


where h is the angular momentum vector in the TEI system 

* TEI TEI 

(computed via eq. 3-15 using R and V from 

SO so 

eqs . 3-1 and 3-2) . ° ° 


The angle from the ascending node to this point, u , is given 

vfi 

by ... 


U CA = tan 


-1 


i h * (ft x r CA ) 


A /\ 

ft • r 


CA 


which can be simplified to 


u CA = tan 


/v A 


-1 ) h • (ft x C) 


7T- 

0 • c 


(3-35) 


where is computed by equation 3-19 using the right ascension 
of the ascending node in the TEI system. 


The TCA, 


tcA' then given by 


II 

< 

u 

M-l 

U CA 

w 


e ca 

2 tan 

-1 

llr+s ' 

m ca = 

e ca " 

e i 

Bin E ca 

fc CA 

m ca + 

X 



tan 


M 



(3-36) 

(3-37) 

(3-38) 

(3-39) 
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where 


x is given by equation 3-30 
M is given by equation 3-28 
e is the orbit eccentricity 
w is the argument of perigee 

This time would then be used in equations 3-32 and 3-33 to update 
the value of C. Equations 3-35 through 3-39 would be repeated 
until the value of t CA between two successive iterations .is 
less than or equal to the time step. At. After the. TCA has been 
found, the corresponding closest approach angle, y , would be 

x-'.rA. 

computed. 


The unit vector from the center of the Earth to the S/C is given 
by 


r = 


cos £2 cos u — sin £2 sin u cos i 
sin £2 cos u + cos' £2 sin u cos i 
sin u sin i 


(3-40) 


The closest approach angle is given by 

-1 


Y CA = cos 


where 


r CA * C CA 


(3-41) 


r CA is 9 iven by equation 3-40 evaluated at u 
^ CA 

C CA is 9 iven by equations 3-32 and 3-33 evaluated 


at t 


CA* 
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If Y is greater than the outer boundary cone angle, y A 
(sec. 3.1), then no AOS or LOS times are possible for this 
'revolution. In this event, the procedure outlined is repeated 
for the next revolution. The initial estimate of TCA for the 
next revolution would be equal to the TCA for the current re- 
volution plus one orbital period. The process continues until 
either y CA < y A or the TCA exceeds the end time of the search, 

t end* If y ca - y a' estiinates of AOS and LOS times are made 
using the method outlined in section 3. 3. 2. 2. 

3 . 3 . 2 . 2 Determining AOS and LOS Times 

After a feasible TCA has been found, the AOS and LOS times can 
be computed by solving the following equations. 

A A 

r • C = cos y A (3-42) 

Substituting equations 3-40 and 3-32 into this expression and 
simplifying yields 

cos (Q - a) cos <j> cos u - sin (fi - a) cos ' d> cos i sin u 

o o 

+ sin sin i sin u = cos y a ( 3— 43 ) | 

O A 

Equation 3-43 cannot be explicitly solved for time for noncircular/ 
nonequatorial orbits. However, this equation can be solved 
numerically using the method of successive substitution. For 
convenience, equation 3-43 can be rewritten as 



A 1 

cos 

u + B 1 

sin u - 

C 1 

(3-44) 

A^ = cos 

*o 

cos 

(fl - a) 



(3-45) 

- sin 

4> 

Y o 

sin 

i - cos 

d> cos 
v o 

iusin (Q - a) 

(3-46) 

C^ = cos 

y a 





(3-47) 
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Assuming a is constant between iterations, equation 3-44 has 
the following solution 


_1 1 

Bl | 

1 , ___-i 1 

c i 

U AOS/LOS 

a i! 

i ■ i 

\k 2 + b i 2 ’ 


The minus sign will produce AOS solutions and the plus sign will 
produce LOS solutions. The procedure is as follows: 


a. Solving for AOS times, t anc 

AUb 

(1) Use initial estimate of t AQS - t CA< 

(2) Compute A^, B^, and using equations 3-33, 3-45, 

3-46, and 3-47. 

(3) Compute u AQS using equation 3-48 with the minus sign. 

(4) Compute revised estimate of AOS time using equations 
3-24 through 3-30. 

(5) If the change in time from the previous iteration is 
less than At, then a solution has been found. Other- 
wise, repeat steps (2) through (5) with the revised time 
estimate. 

b. Solving for LOS time, t T „ 

iiOb 

The solution for LOS time is similar to the procedure for 

AOS time except that 


(1) The initial estimate of LOS time in step a(l) is given 
by 


( 2 ) 

(3) 


^LOS fc CA + (t CA " t AOS ) 


(3-49) 


^LOS is COItl P uted in ste P a (3) by using equation 3-48 
with the plus sign. 

The revised estimate of LOS time is produced in step a (4). 
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The computations discussed in this section are performed for 
each revolution in which y ca was found to be less than or equal 
to y a - In the special case where y CA is exactly equal to Y A > 
the AOS and LOS times are identical and equal to the TCA. Hence 
computation of AOS and LOS times can be bypassed. 
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3.4 ASSUMPTIONS AND LIMITATIONS 

The following assumptions are implicit in the equations presented 
in sections 3 . 1 through 3.3. 

a. Unperturbed orbital motion over the time period t end > t > 
t start. 

b. The effects of polar nutation and precession from t to 

t en d can ke neglected. 
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4.0 FUNCTIONAL OVERVIEW OF AOS AND LOS TIME COMPUTATIONS 

This section provides a functional overview of the ATSVP. The 
overall procedure for determining precise AOS and LOS times is 
presented and discussed. The basic procedure consists of 
predicting the AOS and LOS times using the algorithm presented 
in section 3 and then refining these predictions by performing 
a sequential time search using the equations presented in section 2 . 

The AOS times are defined to be the time point at which the S/C 
first enters the area target or space volume. If the S/C lies 
within the area target or space volume at the beginning of the 
user-specified search period (t s t ar t^ , the first AOS time will 
be set equal to the start time. LOS times are defined to be 
the last time point prior to the S/C exiting the area target 
or space volume. If the S/C lies within the area target or space 
volume at the end of the user-specified search period (t 3 ) t then 
the last LOS time will be set equal to the end time. It should 
be noted that the time increment used to perform the sequential 
time search will limit the accuracy and resolution of the AOS 
and LOS times (e.g., if the time increment is 1 minute, this 
implies that the AOS and LOS times will be determined to the 
nearest minute and that visibility periods of less than 1 minute 
may be skipped) . 

The area targets and space volumes defined in section 1 fall into 
two general categories: 

a. Earth-referenced which: include 

(1) Earth-referenced circles. 

(2) Earth-referenced polygons. 

(3) Earth-referenced space volumes. 

b. Celestial-fixed which include 

( 1 ) celestial-fixed circles. 

( 2 ) celestial-fixed polygons. 

( 3 ) celestial-fixed space volumes. 


4-1 



The logic for computing the AOS and LOS times for each of these 
categories is presented in the following subsections. 


4.1 EARTH-REFERENCED AREA TARGETS ^ND SPACE VOLUMES 


Figure 4-1 provides a functional flowchart for the Earth-refer- 
enced area targets and space volumes. The required inputs are 


a. 

b. 

c . 

d. 

e. 


Target table (either the ground target, ground polygon, or 
ground-fixed space volume table) . 

Target ID. 

Start and end times and time increment (t 


start' 


t , , and 
end 


At, respectively) . 

S/C ephemeris and ephemeris ID. 

RNP matrix and associated epoch time. 


A brief description of the process is provided herein. The 

heading numbers correspond to the numbered blocks on figure 4-1. 

1. The geodetic parameters defining the area target or space 
volume are obtained from the appropriate target table based 
upon the input target ID. Sections 2.1.2, 2.3.2, and 2.5.2 
describe the specific parameters which are required. 

2. The RNP matrix and its associated epoch time are obtained. 

3. The geodetic coordinates of the area target or space volume 
are transformed to the rotating geocentric coordinate system 
(secs. 2.1.2, 2.3.2, and 2.5.2). 

4. The number of visibility periods between t , . and t 

start end 

are determined along with predictions of AOS and LOS times. 
Section 3 describes this process. 

5. A test is performed to determine whether any visibility 
periods were found. If this test is failed, then no AOS/LOS 
times are possible and processing is terminated. 

6. A loop is established which will process each visibility 
period found in step 4. 
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7 . 


The current time, t, is initialized to the AOS time predict- 
ion, t, rtC . . The initial value of the precise AOS time, AOS-, , 
is set to K zero. The AOS/LOS counter, i, is' initialized 
to one. 

A loop is established which will terminate when t exceeds 

the LOS time prediction, t LOg * 

K 

The S/C position vector {at time t in the M50 system) is 
obtained from the ephemeris file based upon the ephemeris ID. 
The S/C position vector is transformed to the rotating geo- 
centric system and computations are performed to determine 
whether the S/C lies within the area target or space volume 
(secs. 2.1.2, 2.3.2, 2.5.2, and app. A). 

Based upon the results of these tests, the S/C visibility 
parameter, V, is set 

V > 0 if the S/C lies within the area target or space 
volume 

V < 0 if the S/C lies exterior to the area target or 
space volume 

11. A test is performed on V. 

If V > 0 

11.1 a further test is performed on AOS.. : to determine 
whether the S/C was also visible during one or more 
of the previous time steps. 

If AOS i > 0 

11,1,1 then the S/C was visible - during one or more 

of the previous time steps. Hence, no transi- 


8 . 

9, 

10 . 


"^Depending upon the target geometry and the S/C groundtrack, 
multiple AOS/LOS times could occur during a single visibility period. 
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tion has occurred. The current time is incre- 
mented by At and the search continues . 

If AOS. < 0 
1 — 

11.1.2 then a transition into the area target or 

space volume has occurred. The current AOS 
time, AOS^, is set equal to the current time. 

The. current time is incremented by At and 
the search continues. 

If V < 0 

11.2 a further test is performed on AOS^ to determine whether 

the S/C was visible during the previous time step. 

If AOS. > 0 
x 

11.2.1 then the S/C was visible during the previous 
time step and a transition out of the area 
target or space volume has occurred. The 
current LOS time, LOS^ , is set equal to the 
time of the previous time step (t - At) . The 
AOS/LOS counter, i, is incremented by one. 

The next AOS time is initialized to zero. The 
current time is incremented by At and the search 
continues. 

11.2.2 the S/C was not visible during the preceding 
time step. Thus no transition has occurred. 

The current time is incremented by At and the 
search continues. 

12. At the completion of the sequential time search, a test is 

made to determine whether the S/C was visible at the end time. 

If the test is true, the last LOS time is set equal to t . 

los k 

4.2 CELESTIAL-FIXED AEEA TARGETS AND SPACE VOLUMES 

Figure 4-2 illustrates the logic flow for the celestial- fixed 
area targets and space volumes.. This. logic is similar to the 
approach presented in section 4.1. The "required inputs are 
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a. Target table (either celestial circles, celestial polygons, 
or celestial-fixed space volumes table) . 

b. Target ID. 

c. Start and end times and time increment t e nd' and 

At, respectively) . 

d. S/C ephemeris and ephemeris ID. 

A brief description of the process is provided herein. The 
heading numbers correspond to the numbered blocks on figure 4-2. 

1. The parameters defining the celestial-fixed area target or 
space volume are obtained from the appropriate target table 
based upon the input target ID. Sections 2.2.2, 2.4.2, and 
2.6.2 describe the specific parameters which are required. 

2. The number of visibility periods between tg-^art and fc end 
are determined along with predictions of AOS and. LOS times. 
Section 3 describes this process. 

3. A test is performed to determine whether any visibility periods 
were found. If- this test is failed, then no AOS/LOS times 

are possible and processing is terminated. 

4. A loop is established which will process each visibility 
period found in step 2. 

5. The current time, t, is initialized to the AOS time prediction, 
t„^„ . The initial value of the precise AOS time, AOS. , is 

aos k A 1 

set to zero. The AOS/LOS counter, i, is initialized to one. 

6. A loop is established which will terminate when t exceeds 
the LOS time prediction, ^lqs * 

7. The S/C position vector (at time t in the M50 system) is 
obtained from the ephemeris file based upon ephemeris ID. 

8. Computations are performed to determine whether the S/C lies 
within the area target or space volume (secs. 2.2.2, 2.4.2, 


Depending upon the target geometry and the S/C ground track, 
multiple” AOS/LOS times could occur during a single visibility 
' period. 
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2.6.2, and app. A). Based upon the results of these tests, 
the S/C visibility parameter, V, is set 

V > 0 if the S/C lies within the area target or 

space volume 

V < 0 if the S/C lies exterior to the area target 

or space volume 

9. A test is performed on V. 

If V > 0 

9.1 a further test is performed on AOS . to determine whether 
the S/C was also visible during one or more of the previous 
time steps . 

If AOS i > 0 

9.1.1 then the S/C was visible during one or more of 
the previous time steps. Hence, no transition 
has occurred. The current time is incremented 
by At and the search continues. 

If AOS. < 0 

l — 

9.1.2 then a transition into the area target- or space 
volume has occurred . The current AOS time , AOS^ , 
is set equal to the current time. The current 

time is incremented by At and the search continues. 

If V < 0 

9.2 a further test ' is performed on AOS^ to determine whether 
the S/C was visible during the previous step. 

If AOS i > 0 

9.2.1 then the S/C was visible during the previous time 
step and a transition out of the area target or 
space volume has occurred. The current LOS time, 

LOS. , is set equal to the time of the previous 
time step (.t - At). The AOS/LOS counter, i, is 
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incremented by one. The next AOS time is initial- 
ized to zero. The current time is incremented by 
At and the search continues. 

9.2.2 the S/C was not visible during the preceding time 
step. Thus, no transition has occurred. The 
current time is incremented by At and the search 
continues . 

10. At the completion of the sequential time search, a test is 
made to determine whether the S/C was visible at the end 
time. If the test is true, the last LOS time is set equal 
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5.0 DETAILED LOGIC FLOW 


Figure 5-1 presents the detailed logic flow for the area targets 
and space volumes processor. This flowchart integrates the 
functional overview presented in section 4 with the equations and 
approach presented in sections 2 and 3 and the appendixes. 

Annotations are provided on the flowchart to describe the comput- 
ations and tests that are performed. In order to simplify 
figure 5-1, several notations have been used to represent corre- 
sponding FORTRAN functions, e.g. , 


4> (I) 
R (I) 


N (J,I) 


K - 



L = 


1 

1 


represents an array of longitudes indexed by, 
the variable I. 

represents a two-dimensional array [i.e., R(3,5)] in 

which the first dimension represents the Cartesian 
components of the vector and the second dimension is 
indexed by the variable I. Hence, operations involving 
vector quantities imply a "DO loop" on the innermost 
dimension. 

represents a three-dimensional array [i.e., N(3,5,3)] 

in which the first dimension represents the Cartesian 
components of the vector and the second and third 
dimensions are indexed by variables J and • I, respectively. 
Operations involving these vector quantities also imply 
a DO loop on the innermost dimension, 
implies summation over the indicated range using the 
variable specified by the lower bound (e.g., the variable 
L will be used to perform the summation over the range 
L = 1 to L = K - 1) . 
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ENTER 



' Obtain target parameters 

a. Earth-referenced circles 
*(1>» X (1) , h(l) , r c 

b. Celestial circles 

a (1) , 6 (1) , y A 

c. Earth-referenced polygon 
NSIDES, NSEG, VO (J) 

(J = 1,2, ...NSEG) 

(I) , A (I), h (I) 

(I = 1,2,.. .NSIDES) 

d. Celestial polygon 
NSIDES, NSEG, VO (J) 

(J = 1,2,. ..NSEG) 

a (I) , 6(1) 

(I = 1,2,. . .NSIDES) 

e. Earth-referenced space 
volume 

NSIDES, NSEG, VO (J) 

(J = 1,2,. . .NSEG) 

(MI), HD, h x , h 2 

(I = 1,2,... NSIDES) 

f. Celestial space volume 
NSIDES, NSEG, VO (J) 

(J = 1,2,. ..NSEG) 
a (I), 6(1), hj., h 2 

(I = 1,2,. . .NSIDES) 


Target 

table 


If \ 
Target is 
E ar th- referenced 


Figure 5-1.- Detailed flowchart. 
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Earth-referenced targets 


Obtain RNP matrix 
and epoch time 

[RNP]™? 


RNP 

matrix 


DO 

I = 1,NSIDES 


If 

Target is a 
circle 


R em 

cs zz . 

F 

\J 

i y 

'cos 2 [4> (I)] + (1 - F) 2 sin 2 [<HD] 

• 


[h (I) + a p ] cos [X (I) ] cos [$(1)] 

r 

Ti 

H 

,[h (I) + a p ] sin [X(D] cos [<£(1)] 



[h (I) + (1 - F) 2 a p ] sin U(I>] 

1 




' |t(l')| 


Compute geocentric 
vectors to each 
vertex 


Compute angular 
radius (in radians) 


Figure 5-1.- Continued. 



Celestial-fixed targets 



Compute inertial 
vectors to each 
vertex 


angular 
to radians 


Figure 5-1.- Continued. 









* $( 1 ) 
.D 


Earth- referenced or 
celestial circles 


CT. 


B 


NSIDES 

■E 

j = -l 


Z $(j) 


NSIDES 


NSIDES- 


Z 


Earth-referenced or 
celestial polygon 


$ = 

&T = 
U1 B 


J = 1 
NSIDES 



J 


Earth-referenced 
space volume 


-> 

C = 


'NSIDES' 

E 

j — i 

NSIDES 




+ h l> 



Celestial 
space volume 


Compute centroid 
of entire area 
target or space 
volume 


Figure 5-1.- Continued 










Compute 
outer 
boundary 
cone angle 
for 

noncircular 

targets 


Figure 5-1.- Continued 






U1 

I 

VO 
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Figure 5-1.- Continued. 
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SUBDIVIDING CONCAVE POLYGONS 


This appendix discusses the procedure for subdividing complex 
concave polygons into two or more simpler convex segments. The 
purpose of this subdivision process is to permit the equations in 
section 2 to be used on a segment-by- segment basis to test for 
containment. The criterion for the S/C (or S/C subsatellite 
point) to be contained in the concave area target or space volume 
is that it must be contained in any one of its segments. For 
convenience, the equations and figures presented in this appendix 
will use a pentagon as an example. 2 However, this approach is 
easily extended to any n-sided polygon. 

Figure A-l presents three examples of concave pentagons. Figure 
A-l (a) illustrates a pentagon having one concave vertex. Figures 
A-l(b) and A-l(c) illustrate pentagons having two concave vertices. 
These pentagons can always be subdivided into triangles by select- 
ing an "appropriate" vertex and connecting nonadjacent vertices 
(fig. A-2) . The maximum number of triangles necessary to completely 
subdivide any arbitrarily shaped polygon is 


N 


seg. 


- n - 2 


max 


(A-l) 


where 


n = number of sides 


A concave polygon is defined to be a polygon that has one or 
more interior vertex angles exceeding 180 degrees. A convex 
polygon is defined to be a polygon that has all of its interior 
angles less than 180 degrees. 

2 

This corresponds to the maximum number of sides specifically 
addressed in the requirements defined in reference 1. 


A-l 



(a) One concave vertex. 


5 



'1 


(c) Two nonadjacent concave vertices. 


Figure A-l.~ Examples of concave pentagons 




Furthermore, the maximum number of interior angles exceeding 180 
degrees can also be determined by noting that the sum of the 
interior angles of the polygon must be equal to the sum of the 
interior angles of all triangles into which it can be subdivided. 
Thus, the sum of the interior vertex angles for any arbitrarily 
shaped polygon is given by 

n 

2 Y ± = (n - 2) 180 (A- 2) 

i=l 

where 

- interior vertex angles of the polygon 

Thus, the maximum number of interior vertex angles exceeding 180 
degrees, y*, is given by 

N(Y*) max = (n - 3) (A-3) 

This equation limits the maximum number of concave vertices for a 
pentagon to two. Figures A-l (b) and A-l(c) illustrate two ex- 
amples. In figure A-l{b), the two concave vertices are adjacent 
to each other. In figure A-l(c), the two concave vertices are non 1 
ad j acent . 

The selection of the "appropriate" vertex to begin the subdivision 
process is highly dependent upon the shape of the polygon and the 
number and relationship of the concave vertices. Also, it is not 
always necessary to subdivide the polygon into triangles. Figure 
A-3 illustrates another method for subdividing the pentagon of 
figure A-l (a) . In this case, the concave pentagon is subdivided 
into a four- sided convex polygon and one triangle. Furthermore, 
figure A-l by no means exhausts' all of the potential pentagon 
shapes that- could be constructed. 


A- 4 




A— 5 


Since the shape of the area targets and space volumes will remain 
static during a mission, it is recommended that the subdivision 
process be performed manually. There are two distinct advantages 
to this approach: 

a. It eliminates the coding and execution of complex subdivision 
logic. 

b. It can be performed once for each concave area target and 
space volume and does not have to be repeated each time AOS 
and LOS times are desired. 


The treatment of concave polygons will place additional require- 
ments on the target tables other than those specifically mentioned 
in reference 1. In addition to the number of sides and coordinates 
for each vertex, the target tables must also contain the following 
for each polygon-shaped target 


N 

seg 


v°i 


- number of segments into which the target is subdivided 
(3 N se g ii 1 for polygons having five or less sides) 

- integers defining the counterclockwise ordering of the 

vertices for each segment (i = l,2,...N g ) 


The use of these additional parameters can best be illustrated by 
example. For figure A- 3, this pentagon is subdivided into two 
segments. The first segment is a four-sided polygon defined by 
vertices 1, 2, 3, and 5. The second segment is a triangle defined 
by vertices 3, 4, and 5. The corresponding parameters for this 
pentagon would be 


N 

seg 


VO 


1 


VO 


2 


2 

1235 (or 2351 or 3512 or 5123) 
345 (or 453 or 534) 


on 


This can easily be performed by plotting the vertex points 
a Mercator projection. 
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Similarly, for figures A-2(b) and A-2(c): 
a. Figure A-2{b) 


b. Figure A-2(c) 


N 

seg 

— 

3 





r*4 

o 

> 

= 

125 

{or 

251 

or 

512) 

CM 

> 

s 

235 

(or 

352 

or 

523) 

V°3 

= 

345 

(or 

453 

or 

534) 

N 

seg 

= 

3 





V °1 

= 

125 

(or 

251 

or 

512) 

V° 2 

= 

245 

(or 

452 

or 

524) 

vo 3 

= 

234 

(or 

342 

or 

423) 


For consistency, this approach can also be used for convex polygons. 
In this case, ^seg would be one and VO.^ would be set to the 

counterclockwise vertex sequence. 


For computational purposes , the number of sides for each segment, 
n^, can be extracted from the vertex ordering integer, VO^, as 

follows 


n^ = highest values of n^ 


where 


where TRUNC 


VO. 

1 


10 


n i~ 


1 


> 0 


(A- 4) 


TRUNC implies integer truncation. 

The vertex numbers, V . , corresponding to each vertex of the i th 

subpolygon can also be extracted from the vertex ordering integer 
as follows 
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n . -1 

^ x 


j-1 

>i-E 

£=1 


V, 10 


n^-£ 


10 1 


3 2,3 ,... ru 
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CONIC INTERSECTIONS 


This appendix presents the equations to compute the intersection 
points between two right circular cones . 1 The required input 
parameters are 

A 

- unit vector along the axis of the first cone 
Y-^ - half-cone angle of the first cone 
V 2 - unit vector along the axis of the second cone 
Y 2 ~ half-cone angle of the second cone 

The quantities to be computed are 

A A 

I-L and I^ - unit vectors to the intersection points of the 
two cones 


Figure B-l illustrates the two right circular cones and the inter- 
section points. The coordinate system for the cones is arbitrary 
and the only restriction is that both v x and V 2 must have a 

common origin and must be expressed in the same system. The first 
step is to determine whether the relative orientation of the two 
cones permits any intersections. This can be determined by com- 
puting the angle between the axes of the two cones. 

_ A A 

y = cos {V 1 • V 2 } (B-l) 


These equations can also be used to compute the intersection 
points between a right circular cone and a plane or the intersection 
points between two planes. 


B-l 




Figure B-l.- Intersection of two cones. 
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The cones will not intersect if 

Y > Y 1 + Y 2 [fig. B-2 (a) ] (B-2) 

or 

Y + smaller {y-p Y 2 > < greater (y-^ Y 2 ^ ' [ fl 9* B-2 (b) ] . (B-3) 

Assuming intersections are possible (i.e., neither equation B-2 
nor B-3 is satisfied) , the next step is to compute the intersection 
vectors. Figure B-3 illustrates the geometry. The intersection 
vectors are symetrically located on either side of the arc connect- 
ing vectors and V 2 . * Furthermore , the intersection vectors 

lie on the perimeter of both cones. Thus, all of the sides of the 

a a A 

spherical triangle connecting vectors V^, V 2 / an< ^ 1^ are 
known 


A A 


side 

a 

(connecting 

V 1 

and 

*1 

) = Y, 

side 

b 

(connecting 

A 

V 2 

and 

H 
< H 

) = y 2 

side 

c 

(connecting 

1— 1 
< > 

and 

A 

V 2 

) - Y 


The angle from side 'c ' to side 'a' can be determined using the 
half- angle formula for spherical triangles. The result is 


B = 2 tan 1 j (B-4) 

sin (S - y 2 ) ) 

where 



S = 1/2 (y 1 + Y 2 + Y) (B-6) 
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Figure B-2.- Conditions when cones do not intersect. 
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B-5 


The dual value of K in equation B-5 produces two values of B 
in equation B-4. These values of B are equal in magnitude but 

, , v -i 

opposite m sign and correspond to the two intersection points-. 
The unit vectors to the two intersection points are given by 


where 


/N 

£ lx £ 2x 

V lx" 


f cos B sin 

Z l,2 = 

% lY l 2 Y 

v iy 


+sin B sin 


_ £ lz £ 2z 

v iz_ 


cos y., 

l * " L 


V 


1 => (V. , V n , V. ) 

lx' ly' lz 


= > U 2x' £ 2y' S '2z ) 


V 1 XV 2 


I V 1 X V 2 


£ 1 > ( 5 - lx , l lz ) 

A A 

= J 2 x Vl 


(B-8) 


(B~9 ) 


^For the special case of B = 0, the two intersection vectors 
are coincident. This physically corresponds to the situation when 
the two cones are tangent. 
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